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I. INTRODUCTION 



This paper continues our study of nonzero ground state entropy, S ({G},q) ^ 0, i.e., 
ground state degeneracy (per site) W({G},q) > 1, where So = fee In W 7 in g-state Potts 
antiferromagnets on various lattices and, more generally, families of graphs {G}. There 
is an interesting connection with graph theory here, since the zero-temperature partition 
function of the above-mentioned g-state Potts antiferromagnet (PAF) on a graph G satisfies 
Z(G,q,T = 0)paf = P(G,q), where P(G,q) is the chromatic polynomial expressing the 
number of ways of coloring the vertices of the graph G PJ with q colors such that no two 
adjacent vertices have the same color [f|- [JJJ|. Thus, 

W{[ lim G }, q) = lim P(G, q) 1/n (1.1) 

n— >oo n— *oo 

where n = v(G) is the number of vertices of G. An example of a substance exhibiting 
nonzero ground state entropy is ice [f§||. Just as complex analysis provides deeper insights 
into real analysis in mathematics, the generalization from q G Z+, to q G C yields a deeper 
understanding of the behavior of W({G},q) for physical (positive integral) q. In general, 
W({G},q) is an analytic function in the q plane except along a certain continuous locus 
of points, which we denote B. In the limit as n — > oo, the locus B forms by means of a 
coalescence of a subset of the zeros of P(G, q) (called chromatic zeros of G) |TTJ. 

In a series of papers |TJ|- |^TJ we have calculated and analyzed W({G},q), both for 



physical values of q (via rigorous upper and lower bounds, large-g series calculations, and 
Monte Carlo measurements) and for the generalization to complex values of q. In the present 
work we construct a variety of homeomorphic families of graphs and give exact calculations of 
the chromatic polynomials, W functions, and resultant nonanalytic loci B. We also compare 
these loci with zeros of the chromatic polynomials for graphs with reasonably large number of 
vertices. Our results extend our previous study, with M. Rocek, of open (planar) strip graphs 
of lattices |3D| via homeomorphic expansion, as will be discussed in detail below. The families 



of graphs that we will construct and study depend on a set of (positive integer) parameters, 
including homeomorphic expansion indices ki, k 2 , and k; the number m of longitudinal edges 
along the strip; and a vector £ = (<7i, ...,cr m ) that describes the allocation of the indices k\ 
and k 2 in the homeomorphic expansions of the various longitudinal edges. As we shall show, 
the chromatic polynomial for such a graph is independent of £ and only depends on the 
first two homeomorphic expansion indices k\ and k 2 through their sum, p = k\ + k 2 . Since 
the number of vertices n = v is a linear function of p, k, and m, there are several ways of 
producing the limit n — » oo (L denotes limit): 

L m : m — ► oo with k ,p fixed (1.2) 
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Lp '■ P — > oo with m , k fixed (1-3) 

and 

Lfc : k —>■ oo with p, m fixed (1-4) 

We shall concentrate on the L m limit here, i.e., the limit of an infinitely long, finite- width 
planar strip, since this is the most interesting of the three limits and since it constitutes 
a homeomorphic generalization of the previous study in Ref. ||20|| . It is natural to begin 
with the simplest such strips, of minimal nontrivial width for this homeomorphic expansion, 
and we shall do so; as will be seen, the results already exhibit considerable richness and 
complexity. We shall also briefly discuss the other two limits, L p and L\.. 
For our results below, it will be convenient to define the polynomial 



> j=0 s=0 ~ / 



where 



a = q-l (1.6) 

and P(C n , q) is the chromatic polynomial for the circuit graph with n vertices: 

P(C n ,q) = a n + (-l) n a (1.7) 

The terms "edge" and "bond" will be used synonymously. 

The organization of the paper is as follows. In section 2 we discuss the connections be- 
tween the structure of the chromatic polynomials, their generating functions, and recursion 
relations that these polynomials satisfy. In section 3 we construct and classify the homeomor- 
phic expansions of strip (chain) graphs. Section 4 contains our calculations of the generating 
function and hence chromatic polynomials for these homeomorphic families of strip graphs. 
This section also includes the calculation of the W function for the L m limit. In section 5 
we present explicit results for the L m limit for various p and k values, and in section 6 we 
discuss and explain general properties of the analytic structure of W for this limit. Section 7 
presents analogous results for the L p and limits, together with a comparison of the three 
limits. Some concluding remarks are given in section 8. Useful properties of are listed in 
Appendix 1. 
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II. CONNECTIONS BETWEEN STRUCTURE OF CHROMATIC POLYNOMIAL, 
GENERATING FUNCTION, AND RECURSION RELATION 



A. Generating Function and Chromatic Polynomial 

A general form for the chromatic polynomial of an n-vertex graph G is 

N a 

P(G,q) = Co(q)+J2 c MMQ)) tjn (2-1) 

where Cj(q) and a,j(q) are certain functions of q. Here the aj(q) and Cj^o(q) are independent 
of n, while Co(q) may contain n-dependent terms, such as (—1)™, but does not grow with 



n like (const.) n with \const.\ > 1. As before [16], we define a term ae(q) as "leading" if it 



dominates the n — > oo limit of P(G,q); in particular, if iV a > 2, then it satisfies |a^(g)| > 1 
and |a^(g)| > |a,,(g)| for j ^ £, so that \W\ = \ae\ tj . The locus B occurs where there is 
a nonanalytic change in W as the leading terms ag in eq. ( |2.1| ) changes. Note that the 
Co(q) term may be absent. For some families of graphs, the Cj(q) and a,j(q) are polynomials. 
However, there are also many families of graphs for which Cj(q) and a,j(q) are not polynomial, 
but instead, are algebraic, functions of q; for these families the property that the chromatic 
polynomial is, in fact, a polynomial of q, is not manifest from the expression fl2.1| ). A number 
of families of this latter type were analyzed in Refs. p0|j21[| ; for these we used a generating 
function method to calculate the chromatic polynomials. 

Here we discuss the connection between our generating function method and the form 
(P-ip. Let a graph G\ be constructed by a certain number of successive additions of a 
repeating graphical subunit H to an initial subgraph J, where I may be the same as H. 
This procedure is similar to the method described in Ref. |2(J |22]| . However, graphs G\ 
constructed in this fashion are general recursive families of graphs and need not be strip 
graphs. The graphs typically also depend also on other parameters, which will be denoted 
{k} but will be suppressed in the notation in this section; for the applications in the present 
paper these will be certain homeomorphic expansion indices, as will be described in greater 
detail below. The total number of vertices is given by 

n = v(G m ) = K\m + k,q (2.2) 

where K\ and k depend on {k}. Applying the deletion-contraction theorem [^3j to the graph 
Gi, one obtains a finite set of linear equations, with m-independent coefficients, relating 
P((Gi) m , q), P((Gi) m _i,g), P((G f ) m ,q), and possibly P((Gy) m _i, q), where G f , with / = 
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2, ..,rif, denotes other recursively defined families of graphs. As was the case for the strip 
graphs discussed in Ref. [p0| , the set of equations is linear because the deletion-contraction 
theorem is a linear equation among chromatic polynomials. The number n/ of families of 
graphs involved in the set of linear equations is finite because the repeating subgraph if is a 
finite graph. This set of linear equations can be solved by using generating functions of the 
form 

oo 

T(G f ,q,x) = P((Gf) m ,q)x m , where f = l,2,..,n f (2.3) 

m=0 

where x is an auxiliary variable. We find that the generating functions T(Gf, q, x) are rational 
functions of the form 

with 

N{G f ,q,x) = Y< A fM^ (2-5) 

and 

V(G f ,q,x) = l + J2hMx 3 (2-6) 

where the Afj and bfj are polynomials in q that depend on the family of graphs Gf, and 
the degrees of the numerator and denominator, as polynomials in the auxiliary variable x, 
are denoted 

d M (G f )=deg x (M(G f )) (2.7) 



dv{Gf) = deg x (V(G f )) (2.8) 
Let us write the denominator of the generating function T(Gf, q, x) in factorized form 

dx> 

V(G,,q,x) = Y[(l-\ fJ (q)x) (2.9) 

5=1 

Henceforth, for brevity, we suppress the G /-dependence in the notation. Using the partial 
fraction identity 
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n 



i=i v J ; i=i 
with the identity, for each j, 



(1 — Xjx 



n 



3 ' l<i<dp; i^j 



(A, - Aj)J 



(2.10) 



1 — A,-£ 

J m=0 



m x m 



(2.11) 



in eq. fl2.4p, with (|2.5D, we obtain the result 



P(G m ,q) = J2 A sFn 



(2.12) 



s=0 



where 



dp 



F r = x; (a^- 1 n 



j = l L 



l<i<<ix>; i^j 



(Aj - Aj)J 



(2.13) 



with the restriction that dj^ < dx>, as is true for all families of graphs considered. Equiva- 
lent ly, 



dr> r dtf 

P(G m ,q) = J2 J>* A f ' ' 

3=1 L s =0 



n 

l<i<dx>;«7^i 



(Aj - Aj) 



A' 



(2.14) 



For example, for dj^ = 1 and dx> = 2, as was the case for a number of strip graphs studied 
in Ref. 0, eq. ([HJ) reduces to 



P{G m ,q) 



(4>Ai+4i) xm (A A 2 + A! 



7^A™ + 



" A 2 



(Ai - A 2 ) 1 (A 2 - Ai) 
For the more complicated case dj^ = 2 and d-® = 3, eq. ( |2.14|) reduces to 

p ^ v _ (AqXI + AiAi + jg) (A A^ + ^A 2 + A 2 ) ^ m 
( m ' gj " (A 1 -A 2 )(A 1 -A 3 ) Al + (A 2 -A 1 )(A 2 -A 3 ) + 
(A Xl + AiX 3 + A 2 ] 



(2.15) 



(A 3 -A 1 )(A 3 -A 2 ) 



Ai 



(2.16) 



and so forth for higher values of djy and dx>- 

Finally, we use eq. Q2.2| ) to express eq. ( |2.12| ) in the form fl2.1|): 



(2.17) 
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where 



N a = d v 



(2.18) 



s=0 



s-1 



n 



(A, - A») 



a; 



«o/«i 



(2.19) 



and 



1 V 

Kl 



(2.20) 



(2.21) 



The advantage of the formulas for the chromatic polynomial given in eq. ( 2.14j) or equiva- 
lently eqs. ( |2.17] )-( |2.21| ) is that they show manifestly the role of the Ay's as the terms dj, and 
thus show that the equation for the continuous locus of points B where W is nonanalytic is 
given by the degeneracy of leading terms 



I A, 



(2.22) 



where "leading" or maximal Aj is defined in the same way as the definition given above of the 
leading term a,j. However, the fact that the chromatic polynomial is, indeed, a polynomial 
in q, is not always manifest in the formulas (12.14 ), ( 2.17|) — (|2.21|) since the c/s and Ay's may 
be algebraic, but nonpolynomial, functions of q. 



B. Generating Function and Recursion Relation 



Writing eq. ( |2.4|) as VY = Af and substituting eq. (|2.3|) , we obtain 



dp 



s=l m=0 

which can be rewritten in the form 



j=0 



m=dj} 



dp 



s=0 



m=0 



s=0 



3=0 



(2.23) 



(2.24) 



Equating the coefficients of the same powers of x on either side of this equation yields the 
recursion relations J2i=o^sP(Gj-s, = Aj for j < d v and Ylt=o b sP(G 'j-siQ) — Aj for 
j > dx>. These can be expressed as the single recursion relation 
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min(j,dxi) 

b s P{G 3 _ s ,q) = A 3 (2.25) 

s=0 

where min(x,y) denotes the minimum of x and y [pi]] . This recursion relation holds for 
arbitrary djy and d-p and involves only one family of graphs G. 



III. CONSTRUCTION OF HOMEOMORPHIC EXPANSIONS OF FAMILIES OF 

STRIP GRAPHS 

We now study homeomorphic expansions of families of lattice strip graphs, including the 
limit of infinite length (with finite width). Homeomorphic classes of graphs have been of 
continuing interest in graph theory ; for some recent theorems in a different direction from 
that of the present work, see Ref . |25| . We have previously used the method of homeomorphic 



expansion to generate a large variety of families of graphs with noncompact W boundaries E 



in the n — > oo limit ||21|| . We recall the definition that two graphs G and H are homeomorphic 
to each other, denoted as G ~ H, if one of them, say H, can be obtained from the other, 
G, by successive insertions of degree-2 vertices on bonds of G J7j. Each such insertion 
subdivides an existing edge of G into two, connected by the inserted degree-2 vertex. This 
process is called homeomorphic expansion and its inverse is called homeomorphic reduction, 
i.e. the successive removal of vertices of degree 2 from a graph H. Clearly, homeomorphic 
expansion of a graph always yields another graph. The inverse is not necessarily true; i.e., 
homeomorphic reduction of a graph can produce a multigraph or pseudograph instead of 
a (proper) graph ||. Here, a multigraph is a finite set of vertices and bonds that, like a 
graph, has no bonds that loop around from a given vertex back to itself but, in contrast to 
a (proper) graph, may have more than one bond connecting two vertices. A pseudograph is 
a finite set of vertices and bonds that may have multiple bonds connecting two vertices and 
may also have looping bonds. For example, consider homeomorphic reduction of a circuit 
graph C r : removing one of the vertices (all of which have degree 2), one goes from C r to 
C r _i, and so forth, until one gets to C3. During this sequence of homeomorphic reductions, 
one remains within the category of graphs. However, the next homeomorphic reduction takes 
C3 to C2, which is a multigraph, not a proper graph. Removing one of the two vertices in 
C 2 produces the pseudograph C\ consisting of a single vertex and a bond that goes out and 
loops back to this vertex. Thus homeomorphism is an equivalence relation on pseudographs. 
This complication will not be relevant for our actual calculations of chromatic polynomials 
here, because we shall only use homeomorphic expansions, not reductions, of graphs, and 
the homeomorphic expansion of a proper graph always yields another proper graph. 



7 



To investigate the effects of homeomorphic expansions, it is natural to consider first the 
simplest case, viz., an open (planar) chain of p-gons, i.e., polygons with p vertices. We shall 
picture the chain as having its longitudinal direction oriented horizontally and its transverse 
direction oriented vertically, and being comprised of m p-gons, and thus m + 1 vertical edges 
(see Fig. 1, which will be discussed further below). 




(a) 



(b) 





(O 



(d) 




(e) 



FIG. 1. Illustration of homeomorphic expansions of open chains of p-gons. (a) {Ch)k lt k 2 ,k,m with 
k\ = &2 = 2, whence p = 4, and k = 3, m = 4; (b) {Ch) kl ^ 2 ^^^ m with ki = 3, &2 = 4, whence p = 7, and 
E = E+, m =4 = (+,+,+,+), k = 3, m = 4; (c) same as (b), but with E = (+,-,-,+); (d) {Ch) kl ,k 2 ,x,k,m 
with k\ = 1, k,2 = 2, whence p = 3, and E = E a / t .. m= 4 = (+, — , +, — ), k = 3, m = 4; (e) same as (d), but 
with E = E + . m= 4. 

The property that the chain is open is equivalently stated as the property that it has open 
or free boundary conditions in the longitudinal direction. Extending the notation used in 
T8| , |20[| , we consider an open chain of m p-gons constructed such that a given p-gon intersects 



with the next p-gon in the chain along one of their mutual vertical edges, and such that 
for the j'th p-gon, there are k\ vertices along either the upper or the lower side from the 
vertex of a transverse (vertical) edge to the vertex of the next transverse edge (including 
these two latter vertices) and (ii) /c 2 vertices along the other (=lower or upper, respectively) 
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longitudinal side, from the vertex of the given transverse edge to the vertex of the next 
transverse edge, including these two latter vertices. If k\ ^ k 2 , then for each of the p-gons in 
the chain, one has a choice as to whether to assign the k\ vertices to the upper longitudinal 
side and k 2 vertices to the lower side, or vice versa. To incorporate this information in the 
k\ 7^ k 2 case, we define a parity vector 

£ = ((7i,(72, ...<r m ) (3.1) 

in which the Oj = + if for the j'th p-gon, the assignment is (upper, lower) = (k\,k 2 ) and 
Oj = — if the assignment is (upper, lower) = (k 2 , ki). We shall denote this strip graph as 

(3.2) 

with the understanding that if k\ — k 2 , S is not defined and can be omitted in the notation: 

(Ch) 

ki,ki,m 

(3.3) 

Thus, (C , /i) 2i 2,m represents the strip (chain) of m squares, with each successive pair of squares 
intersecting on a common edge. It will be convenient to define two special £ vectors: 

E + , m = (+,+, ...,+) (3.4) 

and, for the case of alternating signs, 

^ait.,m — {+,—,+,—, ■■■) (3-5) 

where the dimension m of £ is indicated explicitly here but will sometimes be implicit below. 
Obviously, the (Ch)k ly k 2 ,T,,m graphs are invariant under a reflection about the longitudinal 
axis, which in the k\ ^ k 2 case amounts to a simultaneous reversal of the signs of all of the 
ct/s. Let us define this parity operation on S as 

P(E) = — E (3.6) 

Then 

(Ch) k (Ch) (Ch) (3.7) 

The total number of vertices in each p-gon is 

p = h + k 2 (3.8) 
and the total number of vertices in the chain is 
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n = v(Ch klM ,v,m) = ip - 2)m + 2 (3.9) 

The chromatic polynomial for this open chain (Ch) kl>k2t ^ m is 

P((Ch) klM ^ m , q) = q(q - l)D p (q) m (3.10) 

Observe the important properties that n and the chromatic polynomial P{{Ch) ku k 2 ^,mi q) 
(i) depend on k\ and k 2 only through their sum, p = k\ + k 2 \ and (ii) are independent of E. 
These properties were implicit in the compact notation that we used previously in Ref. [|18 



viz., simply (C/i) Pjm , representing any of the specific p-gon chain graphs (Ch)k lt k 2 ,T,,m such 
that k\ + k 2 satisfies the condition (|3.8|) . This provides an example of the fact that there 
is not a 1-1 correspondence between graphs and chromatic polynomials of graphs; different 
graphs may have the same chromatic polynomial. In the present case, all of the different 
graphs (Ch)ki,k2,T:,rn satisfying the condition ( |3.8|) have the same chromatic polynomial, eq. 
Q3.10| ). If ki ^ k 2 , then, taking account of the overall symmetry ( |3.7|) , this number is given 
by 2 m_1 . Of course, if ki = k 2 , then there is only one such graph. 

For ki > 2 and k 2 > 2, we may view the family P((Ch) kl ,k2,Ti,mi q) as being constructed 
by means of homeomorphic expansion of the strip of m squares, (Ch) 2 , 2 ^ m : one performs 
homeomorphic expansions by inserting k\ — 2 degree-2 vertices on the upper or lower edge 
of the first square and k 2 — 2 degree-2 vertices on the other side (lower or upper, respec- 
tively) of this first square, making another homeomorphic expansion on the second square, 
and so forth along the strip. The resultant homeomorphically expanded strip is precisely 
(C/i)/ Clj fc 2j s,m; that is, denoting this homeomorphic expansion (HE) of the longitudinal edges 
as iJEL fcl _ 2 ,fc 2 -2,s, we have 

HEL kl _ 2:k2 _ 2 ^((Ch) 2t2 ^ m ) = (Ch) kl (3.11) 



In the L m limit of eq. ( |1.2|) , i.e., m — > oo with fixed k\ and k 2 , the resultant nonanalytic 
locus B = 0, i.e., W is analytic in the full complex q plane, and is given by 

W([\im(Ch) kl ,k 2 ,x, m },q) = (D p ) 1/{p - 2) where p = h + k 2 (3.12) 

m— »oo 

We next proceed to consider the more complicated case of homeomorphic expansions of 
the transverse edges on the strip of p-gons, including those shared by adjacent p-gons and 
the two edges on the transverse ends of the strip of squares. We let k denote the number 
of vertices on each of these m + 1 transverse edges, including the pair of vertices belonging 
to the original strip before the homeomorphic expansion. The resultant family of graphs 
obtained by this process of homeomorphic expansion of transverse edges is 
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HET k . 2 ((Ch) ) = (Ch) (3.13) 

An illustration is given in Fig. |I|(a) for the case m = 4 and k = 3. For k\ > 2 and k 2 > 2, one 
may equivalently regard the family (Ch) kl fa^k,m as being constructed by the simultaneous 
combined homeomorphic expansions of the longitudinal and transverse edges of a strip of m 
squares: 

(Ch) HET k _ 2 HEL k ^ 2 , k2 ^((Ch) 2 ,2,m) (3.14) 

Again, the graph is the same if one reflects it about the longitudinal axis: 

(C^)fci,fc2,E,fc,m = (C'^)fc 1 ,fc 2 ,P(E),A;,m = (Ch) fafa ,£,fc, m (3.15) 

One illustration is given in Fig. |I](b) for the case m = 4, k = 3, k\ = 3, and k 2 = 4, whence 
p = 7, and E = S + m=4 = (+, +, +, +). A second illustration is Fig. |l](c) for the same values 
of m, k, k\, and k 2 but S = (+, — , — , +). One can also relate certain of these homeomorphic 



expansions to strip graphs studied in Ref. [20 



(Ch) 2 ^ 2 fa m = H ET k ^ 2 (G sq{L y =2),rn) (3.16) 



(C7l)l,2,E oW .,Jfc,m = H ET k _ 2 (G t (Ly=2), m ) (3-17) 

and 



(Ch) 3>3 fa m = HET k _ 2 (G hc ( Ly=2 )^ m ) (3.18) 

Here, the abbreviations sq, t, and he stand for strips of the square, triangular, and honeycomb 
lattice, and L y denotes the number of vertices in the vertical direction. 

It is also possible for p to be less than 4, in which case the graph does not represent the 
homeomorphic expansion of an original strip of squares. Two illustrations of this are given 
in Fig. [I], both with p = 3. With no loss of generality, we take k\ — 1 for these examples. 
Figs. |l](d,e) show graphs with k\ — 1, k 2 — 2, k — 3, m — 4, and £ = S a ; t . im=4 and S +)m=4 , 
respectively. Note that in the limit m — > oo, the family {Ch)x,2,T,,k,m with £ = (+, +, +) 
contains one vertex of infinite degree. For p = 2 and k — 2, the chain is actually not a 
(proper) graph but a multigraph (see Fig. 0(d)), which we will not consider. For p = 2 and 
k > 3, the chain degenerates according to 

(Ch) p= 2fa m = (Cfr)i,i,fc,m = HEC k - 3 (K 2 + if m ) = -f/fc,m (3.19) 
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which we have previously studied in detail |^T| , p6| (where the notation HEC was used to refer 
to the homeomorphic expansion of the bonds connecting the vertices in the K 2 subgraph to 
the vertices in the K m subgraph). 
Two special cases are 



{Ch) Pt k, 



m=0 



n 



and 



(Ch) Pt k,m=l — Cp+2fc-4 

where T n is a tree graph [^6] and C n is the circuit graph. 



(3.20) 



(3.21) 




(a) 




(c) 



(b) 




(d) 




(e) 




(g) 



FIG. 2. Illustrations of homeomorphic classes for (Ch)k 1 .k 2 .z,k.m with m = 4: (a) k\ > 2, k 2 > 2 (whence 



p > 4), eq. ( p.22[ ); (b) ki = 1, k 2 > 2 (whence p > 3), S = E oit ., eq. ( |3.23| ); (c) k x = 1, k 2 > 2 (whence 



p > 3), S = £ + , one member of the set of classes in eq. ( |3.24| ); (d) k\ = 1, k 2 = 1, so p = 2, eq. ( |3.25|) ; (e) 



p > 4, k = 1, eq. (|3.26|) ; (f) m = 1, eq. ( |3.27l ); (g) m = 0, eq. ( |3.28D 



Thus, the families (Ch) kl ^ 2! Y:,k,m include the following homeomorphic classes, as illus- 
trated in Fig. H (here HR denotes homeomorphic reduction, and in > 1 in eqs. ( |3.22|) to 

(H)): 



12 



{Ch)k x fa > Y,,k,m ~ (Ch)2,2,2,m ~ Gsq(L y =2),m ~ H R{G S q(L y =2), 

for k x > 2 , k 2 > 2 , and k > 2 



(3.22) 



{Ch) 1M ,T, alt .,k,m ~ (C^)i,2,E aW .,2,m = G t (L y =2), m ~ H R(G t ( Ly =2),m) for A; 2 > 2 and > 2 

(3.23) 

{(^^^--^((C^i^J for A; 2 >2, k > 2 , and £ ^ £ ».} (3.24) 



(G/i), 



(C7i)i 



for k > 2 



(3.25) 



(C^)fci,fe 2 ,2,fc=l,m — P| Cp+2fc-4 ~ P| C 2 ~ Cl ' [ P) C 2 

m;Ti m-2;Ti 



Ci for p>4 (3.26) 



(Ch)k lt k 2 ,T:,k,m=l ~ Ci 



(3.27) 



(C/i) fcl 



fc 2 ,S,fc,m=0 ~ ^2 



(3.28) 



The families C 



sq(L y ),m 



and Gt 



(L y ),m 



were defined above, and some relevant graph theory 



definitions are given in [p6" |. The notation {...} in eq. ( |3.24j ) indicates that this set actually 
subsumes multiple homeomorphism classes. In eq. ( |3.26| ), the notation f] m . Tl G means 
the chain of m repeated graphs G, each of which intersects the next in a single vertex 
(= Ti), and the notation Ci • G means the chain in which the Ci pseudograph intersects 
with the rest of the chain, G, at its single vertex. In our introductory discussion about 
homeomorphic classes we mentioned how, in general, homeomorphic reductions of graphs 
yield objects that are not (proper) graphs, but instead are multigraphs or pseudographs. 
Specifically, here, H R(G sq (L y =2),m) in eq. ( |3.22|) , illustrated in Fig. ||( multigraph, as 

are HR(G t{Ly=2) , m ) in eq. flp3|) (Fig. |(b)), HR((Ch\ 2 ^ m ) in eq. (P§ (Fig. |(c)), 

flm-2;Ti ^2 



Ciineq. (g]§) (Fig. |(e)) 
n contrast, the homeomorphic reduction 



and (C/0i,i, 2 ,m m eq. (Fig. @(d)), while C x 

and Ci in eq. ( |3.27| ) (Fig. 0(f)) are pseudographs 
(C/i) fcl ,fc 2 ,s,fc,m=o ~ T 2 in eq. ( ggg ) (Fig. g(g)) does yield a proper graph. 
The number of vertices of (C/i)^,^^,*;,™ is 



n = v{{Ch) klM ^ Am ) = (k + p - 4)m + k 



(3.29) 



and the girth (= length of minimum length circuit) is 
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g((Ch) kl ,k 2 ,B,k,m) = 2k + p-4 



(3.30) 



where in both cases, p = k\ + k 2 as given by eq. (|3.8|). From eq. ( E>.30| ), it follows that for 
the nontrivial case m > 1, the chromatic number [ 37] is (independently of fc) 



x((C70 fcl 



2 if p is even 

3 if p is odd 
(x — 2 for the trivial case m = 0.) Hence, if m > 1, 

P((Ch) kuk2i Y;,k,m, q = 2) = if p is odd 



(3.31) 



(3.32) 



As well as its obvious importance for the chromatic polynomial, the value of x also affects 
the properties of the asymptotic function W, as will be evident in our explicit results below. 
We also recall our notation q c for the maximal finite real point at which W is nonanalytic 
MM- 



IV. GENERATING FUNCTION, CHROMATIC POLYNOMIALS, AND W 
FUNCTION FOR HOMEOMORPHIC FAMILIES OF STRIP GRAPHS 

We now calculate a generating function for the chromatic polynomial of the {CK)k\^2 y T,,k,m 
strip graph. Our method uses the deletion-contraction theorem to derive a set of recursion 
relations involving {Ch)^^^^^ and certain other ancillary families of graphs, as in Ref. 
P0|| . We shall initially denote the generating function as r((C7i)fc lt fc 2) E,fc, Q, x), where x is an 
auxiliary expansion variable, and yields the chromatic polynomial P((Ch)k lt k 2 ,T,,k,m, q) as the 
coefficient of x m in the Taylor series expansion of the generating function about x = 0: 

oo 

ki, k 2,Tj,k,ni} 

m=0 

The m = term is just the chromatic polynomial of the end-rung of the chain, 

P(CVfc,£,fe,m=0, q) = P(T k , q) = q(q - l)^ 1 (4.2) 



As before [ 3D | , we find that the generating function is a rational function in the auxiliary 
variable x and, separately, in the variable q. In particular, we find that the degrees of the 
numerator and denominator, as polynomials in x, are 

d N = 1 (4.3) 
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d v = 2 



(4.4) 



i.e., 

T{{Ch) klM ^ k , q, x) = r = — — — —7—2 (4.5) 

T>((Ch) kl>k2j x tk ,q,x) l + b 1 x + b 2 x 2 

Two important properties are that the generating function and chromatic polynomial for the 
homeomorphically expanded graph (Ch) kltk2 ^ ik>m (i) ar e independent of £ and (ii) depend 
on k\ and &2 only through their sum, as given in eq. (|3] 



{T((Ch) kl ,k 2 ,s,k, q, x), P((Ch) kltk2jS}k , m , q)} = fns.(p, k, m) where p = h + k 2 (4.6) 

Henceforth we shall incorporate this property in our notation, as follows: we shall use the 
symbol (Ch) Pikjm to encompass the full class of families {Ch) kl ^ 2 ,T,,k,m such that k\ + k% = p: 



and similarly, 



and 



{Ch) p ^,m = {(Ch) klk2 Y:,k,m} (4-7) 



q)=P((Ch) kMk , m ,q) (U 



T{(Ch) P)kl q, x) = r((C/i) fclife2jSifc , q, x) (4.9) 



Since, as noted above in eq. ( |3.19| ), for p = 2, the chain reduces to HEC k ^{K 2 + K 



m ) i 



which we have already studied in detail in Ref. |21j| , we restrict ourselves to the cases p > 3 

here. Since for k = 2, the elementary result ( |3.10|) holds, we shall concentrate on the case 
k > 3. 

We calculate 

A = P((Ch) pAm=0 , q) = q(q - l)*" 1 (4.10) 

(independent of p) and 

Ai = ~q 2 (q ~ l)(-l) p D k D k ^ (4.11) 

h = -(-l) p D k -(q- l)D fc+4 ,_ 3 (4.12) 

b 2 = (-iy(q-l) p ~ 1 D k D k _ 1 (4.13) 
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Factorizing the denominator of the generating function as in eq. Q2.9Q gives 

V((Ch) Ptk ,q,x) = (l-\ 1 x)(l-\ 2 x) 
where A12 = — (1/2) [61 ± vR] with the discriminant R = b\ — 46 2 , i-e., 



[-l) p D k + (q- l)D k+p „ 3 ± VR 



where 



R = {D k f + \{q - l)D k+p ^f + 2(-l) p (g - l)D k D k+p _ 3 
-^-lnq-iy-'D.D^ 

The degree of R, as a polynomial in q, is 

deg q (R) = 2(k + p - A) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



and the term containing the highest power of q in R is [(q — l)D k +p-?\ 2 ■ The chromatic 
polynomial is thus given by the general formulas, eq. ( |2.15| ) or ( |2.14| )-( p.21| ) with eqs. ( |4.3| ), 



Cl 



At + ApAx 
Ai — A2 



M \-Jfc/(fc+p-4) 



(4.18) 



a,- = Xj, j = 1, 2, and 



In the T7i — > 00 limit, 



c 2 



Ai + A A 2 
A2 — Ai 



t; 



K\ k + p — 4 



for j = 1,2 



W([ lim (C% fc)TO ],g) = (Xma X ) 1/{k+p - 4) 



(4.19) 



(4.20) 



(4.21) 



where X max refers to the A.,- (j = 1 or 2) in eq. ( |4.15|) with maximal absolute magnitude. 
The continuous locus of nonanalyticities B in W / ([lim m ^ 00 (C/7) Pi fc,m], q) is given by the locus 
of solutions of the equation 



|Ai| — IA2 



(4.22) 



as in the general eq. fl2.22|) . 
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V. EXPLICIT EXAMPLES 



Our exact calculations of the respective continuous nonanalytic loci B are shown for 
p = 3 and p = 4, and for each p, several values of k > 3, in Figs. |3|-|| Further results for 
p = 5 and p = 6 are given in Appendix 2. (As we have discussed earlier, with the ordering 
of the limits in eq. (2.13) of Ref. fl6|], W has no discrete isolated nonanalytic points.) 
Because a chromatic polynomial P(G, q) has real (indeed, integer) coefficients, B has the 
basic property of remaining invariant under the replacement q — > q*\ B(q) = B(q*). For lack 
of space we have only included some of our results; plots for other values of k up to k = 10 
for various values of p are given in Ref. [28|. As the number of vertices of a graph n — > oo, 
the nonanalytic locus B forms by the coalescence of a subset of the zeros of the chromatic 
polynomials (= chromatic zeros) of the graph. Accordingly, it is also of interest to plot these 
chromatic zeros for reasonably long finite-length chains in the figures, to assess how closely 
these approach the asymptotic curves and line segments. The reason that slightly smaller 
values of m are used for families with larger values of k and/or p is that, owing to the relation 
( |3.29| ), these lead to comparably large number of total vertices. 
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FIG. 3. Locus B and analytic structure of the function H^([lim m ^ 00 (C/i) Pi/ ! C ^ m ], g) for p = 3 and k = (a) 



3 (b) 4. For comparison, chromatic zeros are shown for m = 22. 
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FIG. 4. As in Fig. | for p = 3 and k = (a) 7 (b) 8. For comparison, chromatic zeros are shown for 



m = 10. 
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FIG. 5. As in Fig. for p = 4 and k = (a) 3 (b) 4. For comparison, chromatic zeros are shown for m 



(a) 12 (b) 10. 
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FIG. 6. As in Fig. ||for p = 4 and k = (a) 7 (b) 8. For comparison, chromatic zeros are shown for 
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VI. GENERAL PROPERTIES OF W AND B 



We note several general features. In the following it is understood that (i) k > 3, since 



B = for k = 2; and (ii) p > 3, since we have already analyzed the p = 2 case in detail [^T| 



1. For all q e *B, |g| is bounded above, i.e., the loci B are compact and do not extend 
infinitely far from the origin, q = 0. 



2. For a given (p, fc), the nonanalytic locus in the L m limit of eq. (1^2) is comprised of 
arcs and, if and only if p is odd, also a line segment on the positive real q axis. These 
arcs and possible line segment do not enclose any regions in the q plane. 

3. For a number of values of p and k, B includes support for Re(q) < 0, although it never 
includes a line segment on, or a complex arc crossing or touching, the negative real 
axis. Typically for a given p, as one increases k beyond k = 3 to moderate values, the 
left-hand endpoints of the left-most complex-conjugate arcs move into the Re(q) < 
half-plane. Since these complex-conjugate arcs with left-hand endpoints extending into 
the Re(q) < half-plane intersect the imaginary axis in the q plane away from the 
origin, the loci B also have the property that Re(q) = does not necessarily imply 
that q = 0. These statements about B also apply to chromatic zeros for finite-length 
strips, as is evident from the figures. 

4. For even p, the number of different components (arcs) on B is 

Ncomp. = k + p — 4 for even p (6.1) 

Since no two of these arcs touch each other, the number of arc endpoints is 

N e = 2N comp , = 2(k + p-4) for even p (6.2) 

If p is even and (i) k is even, then these arcs form (k + p — 4)/2 complex-conjugate 
pairs, none of which crosses the real axis, so 

p even , k even =^ no q c (6.3) 

(ii) k is odd, then there are (k + p — 5)/2 complex-conjugate pairs of arcs and, in 
addition, a self-conjugate arc crossing the positive real axis at the point q c . B does not 
contain any line segment on the real q axis. The point q c monotonically increases as 
k increases through odd values for fixed even p and, separately, as p increases through 
even values for fixed odd k. We have 
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lim q c (p even) = lim q c (k odd) = 2 (6-4) 

k odd — >oo p even—*oo 

For example, we have 

q c (p = 4, fc = 3) = 1.56984... (6.5) 

g c (p = 4,Jfe = 5) = 1.68233... (6.6) 

g c (p = 6,jfe = 3) = 1.59750... (6.7) 

9c (p = 6,Jfe = 5) = 1.70074... (6.8) 

Since x = 2 for even p, from eq. (|3.31|) , it follows that for finite odd k, where there is 
a q c defined, q c < X- 

5. For odd p, the number of different components on B is 

-k 



N =2 

1 ' comp. ^ 
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p — 4 for odd p (6.9) 



where denotes the integral part of x. The locus B always includes a line segment 
on the positive real q axis given by 

qi<q<Q2 = 2 (6.10) 

Since q<i is the maximal finite real point where W is nonanalytic, we have 

q c {p odd; any k) — 2 (6-H) 

Since by eq. Q3.31] ), the chromatic number is x — 3 for odd p, it follows that 

q c = x — 1 for odd p (6-12) 

For a given odd value of k, the left-hand endpoint of the line segment, qi, is a mono- 
tonically increasing function of p and approaches the limit 

lim qi — 2 for any (6.13) 

p odd — >oo 
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so that in this limit the line segment shrinks to zero and disappears. For a given odd 
value of p, as k increases through odd (even) values from 3 (4) toward infinity, the 
value of the lower endpoint increases (decreases) monotonically toward a limiting value, 
qi(k = oo). Some examples of these limiting values are 

qi(p = 3,k = oo) = ^ (6.14) 
3 

q x {p = 5, k = oo) = - (6.15) 

qi (p = 7, k = oo) = 1.589755.. (6.16) 

qi (p = 9, k = oo) = 1.647799.. (6.17) 

In general, qi(p fixed, k = oo) is a monotonically increasing function of p and ap- 
proaches the limit ( |6.13[ ) as p —>■ oo. Some specific values of qi for odd p and finite k 
include 

(6.18) 
(6.19) 
(6.20) 
(6.21) 

If k is even, no two of the components of B touches any other, so that the number of 
endpoints of the arcs and the line segment is given by 

N e = 2N comp . = 2(k + p-A) for odd p and even k (6.22) 

If k is odd, then two of the complex-conjugate arcs have right-hand ends that meet 
the right-hand side of the line segment at q c = 2 to form a multiple point (in the 



9i (P = 


3, k 


= 3) = 


1.29560 


01 (P = 


3, k 


= 4) = 


1.35321 


01 (P = 


5, k 


= 3) = 


1.43016 


01 (P = 


5, k 


= 4) = 


1.56984 



technical terminology of algebraic geometry |3(J) with valency 3 on the algebraic curve 
constituted by B. Hence, with N comp , = k + p — 5 from eq. ( |6.9|) , we have, for the 
number of endpoints on B, 

N e = 2{N comp , - 1) + 3 = 2k + 2p - 9 for odd p and odd k (6.23) 

In both of these cases of even and odd k, the endpoints and the odd-fc multiple point 
are located at zeros of R and hence the branch points of v^R in eq. ( f4.15|) . 
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One can understand these observed properties of B as follows. As before, we restrict our 
consideration to the cases that we have not already studied, namely p > 3 and k > 3. The 
boundedness of \q\ for q G B follows because the degeneracy equation eq. (|4.22f ) does not 



satisfy the necessary and sufficient condition that we gave in (section IV of) Ref. |18| for B 
to extend infinitely far from the origin, q = 0. Now for large \q\ (and p > 3), the degree in q 
of the two highest-power terms in Ax,2 are the same: 

deg q ((q - l)D k+p _ 3 ) = deg q (^R) = k + p-A (6.24) 

where the term in R with the highest-power in q is simply [{q — l)D k+p _ 3 } 2 . Hence, for A 1; 
where y/~R enters with a + sign, these highest-power terms add, so deg q {\\) = k — p — 4, 
while for A2, where enters with a — sign, they cancel, so the resultant leading power 
of q for large |g| is smaller. Since the leading powers in q are different for Ai and A2, it 
follows that the degeneracy equation cannot be satisfied for arbitrarily large \q\ and hence 
B is bounded in the q plane. The reasoning underlying the structure of B as a set of arcs 



and the line segment is the same as we gave in our analysis of strip graphs in Ref. pD |. This 
also explains the fact that the endpoints of the arcs are the branch points of the term ^/R 
occurring in A 1,2- 

The third feature that we have observed is of interest since for many simple families 
of graphs that we have studied |il||T7|j21"f , including all of the infinitely long, finite- width 



strip graphs analyzed in Ref. P0 |, the respective nonanalytic loci B have the property that 
if q G B, then Re(q) > and the further property that if q G B and Re(q) = 0, then 
q = 0. We had previously encountered loci B that included support for Re(q) < and for 
which Re(q) = did not imply that q = when studying families with noncompact £>'s that 
extended infinitely far from the origin q = [£[11] . However, here we observe this phenomenon 
for simple families with compact loci B that are bounded in the q plane. This behavior of 
the asymptotic loci is also reflected in the chromatic zeros for finite-n graphs, as is evident 
from the figures. Let us denote the minimum value of Re(q) for q G B as 

{q R )min = min[Re(q G B)\ (6.25) 

As is evident from the figures, for p = 3, {q R ) m in has the values 0.3522 for k — 3, —0.1766 for 
k = 4, with similar negative values for other k that we have studied, while for p = 4, (q R ) m in 
has the values 0.2047 for k = 3, —0.0813 for k = 4, with similar negative values for other k 
that we have studied. As will be discussed further below, as k — > 00 for fixed p > 3 or as 
p — ► 00 for fixed k > 3, some arcs tend to approach the unit circle \q — 1| = 1, so the motion 
of the left endpoints of the leftmost arcs in B is not expected to be a monotonic function of 
k for fixed p or of p for fixed k. 
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Note that both the families in Ref. [^TJ with noncompact *B's and the families with 
compact *B's including support for Re(q) < in the present work are consistent with our 
earlier conjectures (la) and (lb) in Ref. [^(J. We recall that conjecture (la) stated that if go 
is a chromatic zero of a regular lattice graph G with no global circuits f29|, then Re(qo) > 
and, furthermore, the only chromatic zero with Re(qo) = is go — itself. Conjecture 
(lb) stated that if one considers the n = v{G) — > oo limit of a regular lattice graph G 
with no global circuits, then the points q G B have the property that Re(q) > 0, and 
furthermore, the only point on B with Re(q) = is the point q = itself. It is true that our 
homeomorphic expansions of strip graphs do not contain any global circuits. However, our 
results are still consistent with the conjectures in Ref. |2(J because, as discussed in Ref. [pi 



none of the families analyzed there was a regular lattice graph, and this is similarly true of 
the homeomorphic families that exhibit chromatic zeros and B with Re(q) < in the present 
work, i.e., those with p > 3 and k > 3. This is clear since the definition of a regular lattice 
graph requires all vertices except possibly those on the boundary to have the same degree (= 
coordination number). But all of the families with p > 3 and k > 3 involve homeomorphic 
expansion and hence have vertices with degree 2 as well as vertices with the usual degree for 
the given lattice (and boundary vertices). Hence, these homeomorphic expansions are not 
regular lattice graphs. 

For the other properties in the list above, it is necessary to deal with the even-p and 
odd-p cases separately. The degree of the discriminant R, as a polynomial in g, is 2(k+p — 4) 
from eq. (|4.17|) , and (aside from the cases p = 2 and k = 2) R has only simple zeros, so 
that each of these corresponds to a branch point of VR. If p and k are both even, R has no 
real zeros; i.e., all branch points of \/R occur in complex-conjugate pairs, connected by the 
arcs of B. This yields the results in eqs. ( |3.1|) and (|6.2|) . If p is even but k is odd, then R 
has no real zeros while b\ has a single real zero, so that at this zero, Ai = (1/2)\^R = — X 2 , 
which evidently satisfies the degeneracy equation ([4.22 ). This explains why for this case, p 
even and k odd, B contains a self-conjugate arc that crosses the real axis. The zero of b\ 
is precisely g c , and this has the monotonicity properties discussed above, approaching the 
limit q c = 2 as p — > oo through even values for fixed odd k or k — > oo through odd values for 
fixed even p, as in eq. (|6.4j) . The other arcs on B form complex-conjugate pairs away from 
the real axis. 

For odd p, we can explain the existence of the line segment that occurs on the real q axis 
by noting that in this case, for gel, since the first two terms in Ai 2 are real and nonzero, 
there will be such a portion of B on the real axis if and only if R < so that \/R~ is pure 
imaginary, so that Ai = r + ir', which is equal in magnitude to A2 = r — ir', where r and r' 
denote real quantities. Now eq. ( |4.17| ) shows that the degree in q of R is even, so that R > 
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for large positive and negative q. For odd p, R has two simple zeros at the positive real points 
qi and q2, so that R < for qi < q < q 2 , which is the requisite condition for the line segment 



3jJ . Furthermore, q 2 = 2 for arbitrary k. This is a consequence of the fact that b\ contains 



a factor of (q — 2) by eq. (|9.7|) and b 2 oc D^Dk-i also contains a factor of (q — 2) by eq. (|9.Gp 



since either k or k — 1 is odd. The values of q\ follow from a straighforward analysis of the 
R = equation; for example, the values of q\ for the limit k — > oo are obtained from the 
equation (1 — a) 2 — 4a p_1 = where a = q — 1; since p = 2£ + 1 is odd, this can be factorized 
as [1 — a — 2o ][1 — a + 2c/] = 0. Because p is odd, it follows from the general eqs. ( j3.31|) 
and (|3.32j ) that x((Ch) p od d,k,m) = 3 and hence P((Ch) p dd,k,m,q = 2) = as in eq. (|3.32| ). 
A priori , a zero of a chromatic polynomial P(G, q) at some value q = go is n °t necessarily 
also a zero of the asymptotic function W; this is true if and only if the multiplicity of this 
zero in P is proportional to the number of vertices n = v(G) as n — > oo. This is, indeed, 
true of the zero at q = 2 in eq. (|3.32|) , as follows from our previous demonstration that for 
odd p, both b% and b 2 , and hence also Ai^, contain a factor of q — 2. Hence, 

W{[ lim (Ch) p oddAm l q) cx (q - 2) W^y ( 6 . 26 ) 

m— >oo 

We can obtain an additional result for the case p = 3. In this case, from eqs. fl4.15|) and 
( [4.161 ), ^ follows that each term proportional to x in M{{Ch)^k,q^ x ) an d 'D((Ch) 3t k,q,x) 
contains a factor of Dk, and Ai,2 contains an overall factor of (-D^) 1 / 2 , since 

A 1|2 (p = 3) = ~ [(g - 2)D k ± v^R] (6.27) 

where 

R(p = 3) = D k [(g - 2) 2 D fc + 4(g - lf/Vi] (6.28) 

Hence, for the nontrivial case m > 1, 

P{{Ch)z, k ,m, q) = at the zeros of £> fe (6.29) 

and, indeed, each of these zeros of the chromatic polynomial have multiplicity proportional 
to n, so that 

W{[ lim (C7i) 3jfe , m ], g) cx (D fc )^) (6.30) 

and 



As gets large for fixed p, one can see from the figures that some arcs are located close 
to the unit circle \q — 1| = 1, while for p > 4, some arcs are oriented almost orthogonally to 
this circle, and, for odd p, the line segment on the real axis approaches the limit 
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qi(k = oo) < q < 2 (6.31) 

where specific values of the lower endpoint, q\{k = oo), were given above for various values 
of odd p. 

One can also consider the limit as p gets large for fixed k; here one sees an alternation 
of features such as the presence (absence) of the line segment for odd (even) p. However, as 
noted above (cf. eq. (|6.13|) ), since lim p oc u ^oo<?i = 2, the line segment disappears in this 
limit. Comparing the figures with different p but the same k, one notices that some of the 
arcs tend to approach the unit circle \q — 1| = 1. 

Another property that we observe concerns the zeros of the chromatic polynomials 
(=chromatic zeros) for finite-m graphs {Ch), p ^,m- For most of the complex arcs the chro- 
matic zeros lie close to, but not, in general, exactly on, these asymptotic curves. In contrast, 
for odd-p families, where B contains a line segment on the positive real axis, q\ < q < 2, the 
chromatic zeros lie exactly on this segment, i.e., are real. A priori, one could have envisioned 
a different situation in which some chromatic zeros that have real parts Re(q) roughly in the 
range q\ < Re(q) < 2 also have imaginary parts that approach zero as m — > oo so that they 
move in to the real axis vertically or obliquely from above and below to form (part or all of) 
the line segment. But instead, the chromatic zeros that merge to form the line segment are 
real even for finite strips. One curious feature that we observe is that in certain cases even 
for complex arcs, the chromatic zeros lie very close to, or even (to within machine precision) 
on the asymptotic curves, such as for the northeast and southeast arcs in Fig. |S|(b) for p = 5, 
k = 4 and the middle, vertically oriented arcs in Fig. ^|(a) for p = 5, k = 5, and Fig. |TT|(a,b) 
for p = 6, k = 5, 6, as well as the complex self-conjugate arcs in Fig. |TU|(a) for p = 6, k = 3 



and Fig. |TT](a) for p = 6, k = 5. Earlier we had encountered similar behavior for certain 
families of graphs and proved a general theorem that showed this rigorously for a class 



of graphs that we called "p- wheels" [I7|. As future work, it would be worthwhile to try to 



construct an analogous proof in the present case. 

Finally, we observe that B for lim m _ +oc ,(C/i) p= 4 i A; = 3 im , shown in Fig. ^|(a) is related to 
the locus B that we calculated previously for the infinitely long open strip of the triangular 
lattice of width L y = 3 []20| (see Fig. Fig. 5(a) of Ref. |2l|) by a simple horizontal translation 
by 1 unit in the q plane: 

B(qhm m ^oo(Ch) p=4 , k=3 , m = B{q') Gt{Lv=3) for q' = q + 1 (6.32) 

This is a consequence of the fact that 

T>((Ch) p= 4, k=3 ,q,x)=V(G t(Ly=3h q',x) for q' = q + 1 (6.33) 

so that 
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A ( CTO P=4ifc=3 j(g) = A Gt(Ll/=3) ,i(g')> J = 1, 2 for g' = q + 1 (6.34) 

(This relation does not hold for the respective numerator functions A and Ai, but these do 
not affect B.) 



VII. COMPARISON OF LOCI B IN THE L m , L p , AND L k LIMITS 

As discussed in the introduction, there are three different ways of sending the number of 
vertices of the graph {Ch) pkm to infinity, as specified in the three limits L m , L p , and L k , eqs. 
(|1.2|)-( |i~4"D . We have concentrated on the L m limit because it is the most interesting and is 
a generalization of our previous study of infinitely long, finite-width strip graphs of various 



lattices in Ref. [20]. In this section we present the analogous, and much simpler results for B 



and W in the other two limits. We recall first that for both the L p and L k limit, the trivial 
case m = leads to B = so we take m > 1. Furthermore, in the L p limit, as was true for 
the L m limit, the case k = 2 is trivial (cf. eq. (|3.10| )) and again yields B = 0, so for the 
L p limit we take k > 3. In passing we note that whereas for the L m limit, the value p = 2 
leads to a noncompact B(q), extending infinitely far from the origin, q = 0, which we have 
analyzed in detail previously pl |, for the L k limit, B is the same for p = 2 and p > 2 so we 



include all p > 2 in this case. 

From our general results for P((Ch) Pt k, m ,q), we find that in the L k limit, k — > oo for 
fixed m > 1 and p > 2, and, separately, in the L p limit, p — > oo for fixed m > 1 and k > 3, 

#L fc =£!L p = g: |g-l| = l (7.1) 

Clearly, there is a fundamental topological difference between B for the and L p limits on 
the one hand, and the respective *B's for the L m limit with various values of p > 3 and k > 2 
on the other hand: the B's for the and L p limits divide the complex q plane into two 
regions, whereas the B's for the L k limit with various p and k do not. Another important 
difference is that B in the L k limit is actually independent of the specific values of the other 
parameters, m and p, provided that they satisfy m > 1 and p > 2, and, similarly, B in the 
L p limit is independent of the specific values of m and k, provided that they satisfy m > 1 
and k > 3; in contrast, the locus B in the L m limit does depend on the specific values of p 
and k provided that p > 2 and k > 3. 



Defining, as in our earlier work [16-18,21], the region Ri as the region including the 



positive real axis for large q and extending down to the maximal finite real point, q c , where 
W is nonanalytic, we identify 
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g c = 2 for B Lk , B Lp (7.2) 

and 

R 1 = q: |g-l|>l for B Lk , B Lp (7.3) 

The other region is the complement of Ri, 

R 2 = q: |g-l|<l for B Lk , B Lp (7.4) 

For the nontrivial case m > 1 and the respective limits L k with p > 2 and L p with k > 3, 
we find 

W{{Ch) pAm -L k) q) = W((Ch) pAm ;L p ,q) = q- 1 for q E R x (7.5) 



As we have discussed before fig] , in regions other than R\, one can only determine the 
magnitude of W; we obtain 

\W((Ch) pAm ,L kl q)\ = \W((Ch) pAm ,L p ,q)\ = l for q G R 2 (7.6) 

As simple examples of the cases of (i) fixed k, m and variable p and (ii) fixed p, m and 
variable k, one may take, respectively, (i) k — 3 and m = 1, and (ii) p = 3, m = 1. In these 
cases, the general results (|3.21| ) and ( p..7| ) yield 

P((C/i) Plfe=3>m =i, g) = P(C P+2 , g) = a p+2 + (-l) p a (7.7) 

^((^^=3,^=1, g) = P(C 2k -i, q) = a 2h - 1 + (-l) 2k - l a (7.8) 

where a = q — 1. One sees the results (|7.1Q , (|7.5|), and (|7.6|) immediately. For these cases, 
as we showed previously fll6,P7|, the zeros of the chromatic polynomial lie precisely on the 



asymptotic locus (|7. 1| ) even for finite p. As a more complicated and generic example of 
category (ii), fixed p, m, and variable k, we take p = 3 and m = 2. In Fig. [7| we show the 
zeros of this chromatic polynomial for k = 25, in comparison with the asymptotic locus B 
given by eq. (|7.1| ). One sees that the chromatic zeros lie close to B. 
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FIG. 7. Locus B and analytic structure of the function W([lmik^oc>(Ch)p^,m], <?) for p — 3 and m = 2. 
For comparison, chromatic zeros are shown for k = 25. 

It is interesting to inquire whether, if we first take m —>■ oo with fixed p and k to obtain 
£>L m , and then take, say, p — > oo, the resulting locus is equal or unequal to that which 
we obtain by first taking p — > oo for fixed k and m, thereby getting Bl p , and then taking 
m — > oo. Indeed, one can formulate three such questions: 



1. Is linifc^oo Bg ] L m equal or unequal to linim^oo Bq- L k for fixed p ? 

2. Is linip^oo Bg- L m equal or unequal to linim^oo Bq- l p for fixed k ? 

3. Is linip^oo Bq- L k equal or unequal to lirn^oo Bg-, l p for fixed m ? 

As discussed above, we shall restrict ourselves here to the nontrivial cases m > 1 and 
k > 3, and, since we have already discussed the noncompact case where p = 2 for the L m 
limit, we also restrict to p > 3 for this L m limit. We first observe that from eq. ( |7. 1|) , 
it follows immediately that the two limits in item (3) above do commute. As regards the 
k — > oo limit of £>L m , our results in the previous section show that for odd p, Bi m includes 



the line segment (|6.31|) in this limit, and therefore the limits in item (1) above are definitely 
not equal. Concerning the limits in (2) above, although some of the arcs forming B in the 
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L m limit approach the unit circle ( [7.1|) as p — > oo, others appear to be roughly orthogonal 
to this circle, so the limits in (2) could only be equal if these roughly orthogonally oriented 
arcs shrink and disappear as p — > oo. Thus, the families of graphs studied here exhibit 
both noncommuting limits of nonanalytic loci as certain parameters are taken to infinity in 
different orders, and an example of commuting limits (indeed, in this case, the two loci are 
identical, as noted in eq. (|7. 1| ) . 

VIII. CONCLUSIONS 

In this paper, using a generating function method, we have presented exact calcula- 
tions of the chromatic polynomials for a variety of homeomorphic families of strip graphs, 
(C/i)fc li fc 2i s,fc,m- We have studied each of the ways of sending the total number of vertices 
n — > oo and have determined the respective loci B where W is nonanalytic. In particular, 
we have studied the case of infinitely long strips, i.e., m — > oo, with fixed fci, /c 2 , X (of which 
only the combination p = k\ + hi affects £>), and k. Our results exhibit considerable richness 
and complexity in the forms of B. In the future, it would be of interest to continue the 
calculations of chromatic polynomials, W functions, and the the loci B for other types of 
homeomorphic families of graphs. Our findings yield further insights into the behavior of 
the nonzero ground state entropy of g-state Potts antiferromagnets on various graphs. 

This research was supported in part by the NSF grant PHY-97-22101. 
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IX. APPENDIX 1 



In this Appendix we gather together some convenient formulas concerning the function 
D k , defined in eq. ( |1.5|) in the text. From the definition, it follows that k > 1. Some special 
values are D\ = 0, D 2 = 1, and D% = q — 2. The deg polynomial in g, of D k is 

clearly 

deg q {D k ) — k — 2 for k>2 (9.1) 

In general, one has 

£>*(? = 0) = (-1)*(A; - 1) (9.2) 

and 

D k (q = 1) = (-l) fc for fc > 2 (9.3) 

Since for k even, the circuit graph C k is bipartite, which is equivalent to the fact that the 
chromatic number x(Cfc even) = 2, and thus P{C k even , q = 2) = 2, it follows that 

£fc even{q = 2) = 1 (9.4) 

Since for k odd, x(Cfc od d) = 3 and P(C k odd , q = 2) = 0, we have 

D k odd {q = 2) = (9.5) 

This zero results from a linear factor, i.e., 

Dk odd = {q- 2)pol(q) (9.6) 

where pol(q) is a polynomial of degree k — 3 in q with po/(g = 2) 7^ 0. Further, for p > 3 
and k > 2, 

D k - aD k+2n = (q- 2)pol'{q) (9.7) 

where pol'(q) is another polynomial. Some identities that we have derived and used for our 
calculations are listed below: 

D k - aD k ^ = {-if (9.8) 
or equivalently, D k+ i — D k = D 3 D k + (— l) k+1 , and 
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D 2fc _! = [qD k - 2(-l) k }D k 



(9.9) 



We note some further important properties. If k > 4 is even, then the derivative dD k /dq 
vanishes at a single real value of q, which we denote q m , and d 2 D k / dq 2 \ q=qm > 0, so that this 
is an absolute mininum of D k even (the subscript m on q m refers to this minimum). Note 
that the restriction on k is imposed since D2 = 1 so dD2/dq = 0. Furthermore, 



Ac even{q = Qm) > (9.10) 



so that 



If k > 4 is even, then D k > V q G R (9.11) 
We find that q m increases monotonically from 

q m (k = A) = ^ (9.12) 

approaching 

lim q m = 2 (9.13) 

If k > 3 is odd, then ^ > VgGl (9.14) 

dq 

so that in this case D k has no critical points (zeros of dDk/dq). Together with eq. ( |9.5|) , 
this implies that if k is odd and k > 3, then D k has a single real zero, at q — 2. The proofs 
follow immediately from the definition of D^- 



X. APPENDIX 2 

In this appendix we show exact calculations of the respective continuous nonanalytic loci 
B are shown for two higher values of p, viz, p = 5 and p = 6, and for each p, several values 
of k > 3. These results are discussed in the text. 
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FIG. 8. As in Fig. | for p = 5 and k = (a) 3 (b) 4. For comparison, chromatic zeros are shown for 



m = 12. 
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FIG. 9. As in Fig. for p = 5 and k = (a) 5 (b) 6. For comparison, chromatic zeros are shown for m 



(a) 12 (b) 10. 
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FIG. 10. As in Fig. | for p = 6 and k = (a) 3 (b) 4. For comparison, chromatic zeros are shown for 



m = 10. 
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FIG. 11. As in Fig. | for p = 6 and k = (a) 5 (b) 6. For comparison, chromatic zeros are shown for 



(a) 8 (b) 6. 
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